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© Safe Onboard Guidance and Control Under 
Probabilistic Uncertainty 

NASA’s Jet Propulsion Laboratory, Pasadena, California 

ing optimization is convex. Empirical re- 
sults show that the approach is many or- 
ders of magnitude less conservative than 
existing set conversion techniques, for a 
small penalty in computation time. 

This work was done by Lars James Black- 
more of Caltech for NASA’s Jet Propulsion 
Laboratory. For more information, contact 
iaoffice@jpl.nasa.gov. NPO-46155 


An algorithm was developed that de- 
termines the fuel-optimal spacecraft 
guidance trajectory that takes into ac- 
count uncertainty, in order to guarantee 
that mission safety constraints are satis- 
fied with the required probability. The 
algorithm uses convex optimization to 
solve for the optimal trajectory. Convex 
optimization is amenable to onboard so- 


lution due to its excellent convergence 
properties. 

The algorithm is novel because, unlike 
prior approaches, it does not require 
time-consuming evaluation of multivari- 
ate probability densities. Instead, it uses a 
new mathematical bounding approach 
to ensure that probability constraints are 
satisfied, and it is shown that the result- 


© General Tool for Evaluating High-Contrast Coronagraphic 
Telescope Performance Error Budgets 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


The Coronagraph Performance Error 
Budget (CPEB) tool automates many of 
the key steps required to evaluate the 
scattered starlight contrast in the dark 
hole of a space-based coronagraph. The 
tool uses a Code V prescription of the 
optical train, and uses MATLAB pro- 
grams to call ray-trace code that gener- 
ates linear beam-walk and aberration 
sensitivity matrices for motions of the 
optical elements and line-of-sight point- 
ing, with and without controlled fine- 
steering mirrors (FSMs). The sensitivity 
matrices are imported by macros into 
Excel 2007, where the error budget is 
evaluated. The user specifies the partic- 
ular optics of interest, and chooses the 
quality of each optic from a predefined 
set of PSDs. The spreadsheet creates a 
nominal set of thermal and jitter mo- 
tions, and combines that with the sensi- 
tivity matrices to generate an error 
budget for the system. 


CPEB also contains a combination of 
form and ActiveX controls with Visual 
Basic for Applications code to allow for 
user interaction in which the user can per- 
form trade studies such as changing engi- 
neering requirements, and identifying 
and isolating stringent requirements. It 
contains summary tables and graphics 
that can be instantly used for reporting re- 
sults in view graphs. 

The entire process to obtain a coron- 
agraphic telescope performance error 
budget has been automated into three 
stages: conversion of optical prescrip- 
tion from Zemax or Code V to MACOS 
(in-house optical modeling and analysis 
tool), a linear models process, and an 
error budget tool process. The first 
process was improved by developing a 
MATLAB package based on the Class 
Constructor Method with a number of 
user-defined functions that allow the 
user to modify the MACOS optical pre- 


scription. The second process was mod- 
ified by creating a MATLAB package 
that contains user-defined functions 
that automate the process. The user in- 
terfaces with the process by utilizing an 
initialization file where the user defines 
the parameters of the linear model 
computations. Other than this, the 
process is fully automated. The third 
process was developed based on the 
Terrestrial Planet Finder coronagraph 
Error Budget Tool, but was fully auto- 
mated by using VBA code, form, and 
ActiveX controls. 

This work was done by Luis F. Marchen of 
Caltech for NASA’s Jet Propulsion Laboratory. 
For more information, contact iaoffice @ 
jpl.nasa.gov. 

This software is available for commer- 
cial licensing. Please contact Daniel Brod- 
erick of the California Institute of Technol- 
ogy at danielb@caltech.edu. Refer to 
NPO-4 7220. 


© Hidden Statistics of Schrodinger Equation 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Work was carried out in determination 
of the mathematical origin of random- 
ness in quantum mechanics and creating 
a hidden statistics of Schrodinger equa- 
tion; i.e., to expose the transitional sto- 


chastic process as a “bridge” to the quan- 
tum world. The governing equations of 
hidden statistics would preserve such 
properties of quantum physics as super- 
position, entanglement, and direct-prod- 


uct decomposability while allowing one 
to measure its state variables using classi- 
cal methods. In other words, such a sys- 
tem would reinforce the advantages and 
minimize the limitations of both quan- 
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turn and classical aspects, and therefore, 
it will be useful for implementation of 
quantum computing. 

Recent advances in quantum informa- 
tion theory have inspired an explosion 
of interest in new quantum algorithms 
for solving hard computational prob- 
lems. Three basic “non-classical” proper- 
ties of quantum mechanics — superposi- 
tion, entanglement, and direct 
tensor-product decomposability — were 
main reasons for optimism about capa- 
bilities of quantum computers and 


quantum communications as well as for 
a new approach to cryptography. How- 
ever, one major problem is keeping the 
components of a quantum computer in 
a coherent state, as the slightest interac- 
tion with the external world would cause 
the system to decohere. Another prob- 
lem is measurement: by the laws of quan- 
tum mechanics, a measurement yields a 
random and incomplete answer, and it 
destroys the stored state. 

This proposed reinterpretation of 
quantum formalism opens up new ad- 


vantages of quantum computers: if the 
Madelung equations are implemented 
on a classical scale (using, for instance, 
electrical circuits or optical devices) , all 
the quantum effects important for com- 
putations would be preserved; at the 
same time, the problems associated with 
decoherence and measurement would 
be removed. 

This work was done by Michail Zak of Caltech 
for NASA’s Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl.nasa.gov. 
NPO-46731 


© Optimal Padding for the Two-Dimensional Fast Fourier Transform 

Appending data to an optimum length decreases computing runtime. 

Goddard Space Right Center, Greenbelt, Maryland 


One-dimensional Fast Fourier Trans- 
form (FFT) operations work fastest on 
grids whose size is divisible by a power of 
two. Because of this, padding grids (that 
are not already sized to a power of two) 
so that their size is the next highest 
power of two can speed up operations. 
While this works well for one-dimen- 
sional grids, it does not work well for 
two-dimensional grids. 

For a two-dimensional grid, there are 
certain pad sizes that work better than 
others. Therefore, the need exists to 
generalize a strategy for determining 
optimal pad sizes. There are three steps 
in the FFT algorithm. The first is to per- 
form a one-dimensional transform on 
each row in the grid. The second step is 
to transpose the resulting matrix. The 
third step is to perform a one-dimen- 
sional transform on each row in the re- 
sulting grid. Steps one and three both 


benefit from padding the row to the 
next highest power of two, but the sec- 
ond step needs a novel approach. 

An algorithm was developed that struck 
a balance between optimizing the grid 
pad size with prime factors that are small 
(which are optimal for one-dimensional 
operations), and with prime factors that 
are large (which are optimal for two-di- 
mensional operations). This algorithm 
optimizes based on average run times, 
and is not fine-tuned for any specific ap- 
plication. It increases the amount of times 
that processor-requested data is found in 
the set-associative processor cache. Cache 
retrievals are 4-10 times faster than con- 
ventional memory retrievals. 

The tested implementation of the algo- 
rithm resulted in faster execution times 
on all platforms tested, but with varying 
sized grids. This is because various com- 
puter architectures process commands 


differently. The test grid was 512x512. 
Using a 540x540 grid on a Pentium V 
processor, the code ran 30 percent faster. 
On a PowerPC, a 256x256 grid worked 
best. A Core2Duo computer preferred ei- 
ther a 1040x1040 (15 percent faster) or a 
1008x1008 (30 percent faster) grid. 

There are many industries that can 
benefit from this algorithm, including 
optics, image-processing, signal-process- 
ing, and engineering applications. 

This work was done by Bruce H. Dean, 
David L. Aronstein, and Jeffery S. Smith of 
Goddard Space Flight Center. Further infor- 
mation is con tained in a TSP (see page 1 ). 

This invention is owned by NASA, and a 
patent application has been filed. Inquiries 
concerning nonexclusive or exclusive license 
for its commercial development should be ad- 
dressed to the Patent Counsel, Goddard Space 
Flight Center, (301 ) 286-7351. Refer to GSC- 
15678-1. 


© Spatial Query for Planetary Data 

This technology is extensible to Earth science and satellite monitoring and surveillance. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Science investigators need to quickly 
and effectively assess past observations of 
specific locations on a planetary surface. 
This innovation involves a location- 
based search technology that was 
adapted and applied to planetary sci- 
ence data to support a spatial query ca- 
pability for mission operations software. 

Conventional databases of planetary 
datasets are indexed and searchable by 
various metadata, such as acquisition time, 


phase of mission, and target. Searching 
these datasets will produce enormous 
datasets that are difficult, or impractical, 
to browse through to identify observations 
of very specific targets. For queries at spe- 
cific locations, it is fundamentally more ef- 
ficient to specify the location as the target 
of the query; and to have the database 
search based on the location of the data 
rather than metadata that is only indirectly 
or tangentially related to location. 


High-performance location-based 
searching requires the use of spatial data 
structures for database organization. Spa- 
tial data structures are designed to organ- 
ize datasets based on their coordinates in 
a way that is optimized for location-based 
retrieval. The particular spatial data struc- 
ture that was adapted for planetary data 
search is the R+ tree. The R+ tree 
arranges data as a set of nodes that repre- 
sents bounding rectangles. Every leaf 


32 


NASA Tech Briefs, April 201 1 


